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Abstract 

A brief review is provided of the use of the Square Wave Method 
(SWM) in the held of signal and image analysis and it is specihed how 
results thus obtained are expressed using the Square Wave Transform 
(SWT), in the frequency domain. To illustrate the new approach intro¬ 
duced in this held, the results of two cases are analyzed: a) a sequence 
of samples (that is, measured values) of an electromyographic recording; 
and b) the classic image of Lenna. 
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1 Introduction 

It was previously shown how a new method, the Square Wave Method (SWM), 
for the analysis of signals depending on one variable [1] can be presented in the 
frequency domain by using a mathematical tool called Square Wave Transform 
(SWT) [2] [3]. The SWM was then generalized quite naturally and directly for 
image analysis [3]. 

The objectives of this paper are the following: 

1. To provide a brief review of the use of the SWM for the analysis of signals 
and specify the relations existing between a) the sampling frequency /g, 
with which the successive values of recordings of biomedical signals (such 
as those of an electrocardiogram, electromyogram or electroencephalo¬ 
gram) are measured, and b) the frequencies /i, / 2 , • • •, /n? corresponding 
respectively to the different trains of square waves S'!, £' 2 ,..., 5'n obtained 
using the SWM; 
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2. To indicate how it also is possible to present in the frequency domain using 
the SWT, the results of the analysis of images obtained with the SWM. 

The application of the SWM in the field of signal and image analysis is 
exemplified with the results of an analysis of a) a sequence of samples (that is, 
measured values from an electromyographic recording); and b) the classic image 
of Lenna, using the SWT [5] . 

2 Analysis of a Function of One Variable 

Consider a function of time (t), in the interval At, satisfying the conditions of 
Dirichlet [6]: 


/(t) = (6 — t)(2 cos(27r4t) + 5 cos(27r6t)) 0<t<4s (1) 

Suppose that the time interval in which the function characterized in equa¬ 
tion 0 will be analyzed (At = 4 s) has been divided into 18 equal sub-intervals. 
In this case, it will be seen that function ( 1 ) can be approximated in At, using 
the sum of the parts corresponding to At of 18 trains of square waves. These 
trains of square waves will be called Si^ S 2 , , Sis] the “S'” being based on 

the word “square” in the expression “train of square waves”. 

If At has been divided into 100 equal sub-intervals, the approximation to the 
function 0 in interval At will be carried out by adding the parts corresponding 
to At of 100 trains of square waves: Si, S 2 , S 3 ,..., Sioo- In general, if At is 
divided into any natural number n of equal sub-intervals, the approximation 
in At to function 0 will be obtained by adding the parts corresponding to 
At of n trains of square waves: Si, S 2 , S 3 ,..., S^. The Square Wave Method 
(SWM) described in this section makes it possible to determine those trains 
of square waves unambiguously. Therefore, each Si (where i = 1,2, ...,n) of 
those trains of square waves will be characterized by a specific frequency fi 
(i.e., consideration is given to the number of waves in the train of square waves, 
which is contained in the unit of time 1 s) and a particular coefficient Q, whose 
absolute value is the amplitude of the corresponding train. 

The function f{t) specified in 0 is shown in figure 
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Figure 1: f{t) = (6 — t)(2 cos(27r4t) + 5 cos(27r6t)) 0 < t < 4 s. 

For the case considered here, n = 18, a description will be provided below of 
how the frequencies fi (where i = 1,2,...,18) and the values of the coefficients 
Ci (where i = l,2,...,18) corresponding to the different trains of square waves 
Si (where i = 1, 2 ,..., 18) are determined; see figure 
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Figure 2 : How to apply the SWM to the analysis of the function represented in 
figure(See indications in text.) 


The first row of figure(with coefficients Ci) represents half a square wave, 
the first semi-wave of the train of square waves Si. The frequency of Si (i.e., 
fi) is clearly equal to the number of square waves per unit of time (Is). To 
obtain /i, the part of Si which occupies At (the half-wave) is divided by At. 
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To compute / 2 , note that At is occupied by the sum of that half-wave of the 
train of square waves S 2 and the fraction ^ of the second semi-wave of the first 
square wave of S 2 • That fraction is represented by the symbol C 2 in the second 
row of figure Thus the following value is obtained for / 2 : 


/2 = 



2 ^ 17) _ 1 TL = ^ ^ 

At 2 ’ At 2At VI 8 - ly 


18 


18-1 
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To compute /s, note that At is occupied by the sum of that half-wave of the 
train of square waves £'3 and the fraction ^ of the second semi-wave of the first 
square wave of S^. This fraction (in the third row of figure]^ is represented by 
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the sequence of symbols —Cs'—Cs. Therefore, the following value is obtained 
for /a: 
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Observe that any of the 18 values of fi, where i = 1, 2,..., 18, can be com¬ 
puted with the following equation: 


fi = 


1 


18 


18 - (i - 1) 




i = 1,2,...,18 


In general, if the interval At, whose value, of course, may be different from 
4 s, is divided into n equal sub-intervals, the frequencies corresponding to each 
of the n trains of square waves are as follows: 




1 


n \ 


-.-1. 


t = 1, 2,..., n 


( 2 ) 


It has been explained how to compute each fi corresponding to each 5'^, 
where i = 1, 2,..., 18, for the case of the approximation to /(t) specified in Q 
when dividing At into 18 equal sub-intervals (n = 18). Indications will now be 
given on how to compute the Ci for each Si. 

The vertical arrow pointing down at the right of figure 2 indicates how to 
add the terms corresponding to each of the 18 sub-intervals of At. Thus, to 
obtain the values of the coefficients Ci, C 2 ,..., Cu and C 18 , corresponding to 
S'!, 5'2, • • - 5 Si 7 and 5'i8, the following system of linear equations must be solved. 
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In the preceding system of linear algebraic equations hi, V 2 ,..., Vi 7 and 
V 18 are the values for f(t) as specified in 0 at the midpoints of the first, second, 
third, ..., seventeenth and eighteenth sub-intervals, respectively, of interval At, 
in which f(t) is analyzed. It follows that the values Vi (where i = 1,2,3,. ..,17 
and 18) can be computed given that /(t) has been specified in 0. These values 
are as follows; 


Vi = -34.5484836 
V 2 = 30.6666667 
V 3 = -16.0256827 
V 4 = -6.9904692 
V 5 = 49.0000000 
Ve = -6.5602864 
V 7 = -14.1121683 
Vs = 25.3333333 
Vg = -26.7629098 


Vio = -25.7897131 
Vii = 22.6666667 
V 12 = -11.7202754 
Vi3 = -5.0546469 
Vu = 35.0000000 
Vi5 = -4.6244642 
V 16 = -9.8067610 
Vi7 = 17.3333333 
V 18 = -18.0041393 


Each of the 18 values of Vi, where i = 1,2,..., 18, has been computed with 
a precision of seven decimal digits. 

The 18 unknowns of the systems of equations specified in 0 are Ci, C 2 ,..., 
( 7 i 7 , and Cig. Thus \Ci\ refers to the amplitude of the train of square waves Si, 
where i = 1, 2,..., 18. The (constant) value of each positive square semi-wave 
of the train of square waves Si is \Ci\ and the (constant) value of each negative 
square semi-wave of that Si is —\Ci\. 

The system of equations 0 has been solved by using LAPACK [7] , and the 
following results were obtained for the unknowns: 


Cl = 117.12980 
C 2 = 50.27631 
C 3 = -210.98830 
C 4 = -53.27896 
C 5 = 9.35088 
Ce = 12.58025 
C 7 = 61.27212 
Cg = 49.80105 
C 9 = 12.81335 


Cio = 4.03101 
Cii = -85.68506 
C 12 = 12.88482 
Ci3 = 8.51973 
Ci4 = 60.38772 
Ci5 = -69.86421 
C 16 = 28.08997 
Ci7 = -9.26140 
C 18 = -32.60758 


The trains of square waves 81 , 82 , Ss, ■■■, 817 and Sig have been shown for 
interval At in figures |33.1[ | 33 . 2 [ |33.3| ..., |33.18[ respectively. 
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The approximation obtained for f{t) (as specified in (iT), in interval At, by 
adding the 18 trains of square waves) is displayed in figure 3 



Figure 4: The dashed line indicates the approximation to /(t), specified in Q, 
18 

by 


If one requires a better approximation to /(t), by adding the trains of square 
waves, then At should be divided into a larger number of equal sub-intervals. 
The higher the number of sub-intervals, the better the approximation. 

Suppose that interval At is divided into n sub-intervals of equal duration. In 
equation it was specified how to compute each fi (where i = l,2,3,...,n) 
for each train of square waves S'!, £' 2 , S's,..., 5'n, which must be added in At 
to obtain, in that interval, the corresponding approximation to f{t) specified 
in 0 . 

To obtain the coefficients Ci, (^ 2 , C 3 ,..., corresponding respectively to 
those square waves, a system of linear algebraic equations must be solved. This 
system can be obtained by using the same type of approach as that used to 
obtain the system of equations specified in 0 . 

Approximations to f{t) specified in equation § when dividing At into 100 
and into 1000 intervals respectively, are displayed in |(5.1)| and (5.2)| in figure 5. 

The SWM cannot be considered a branch of Fourier analysis; the trains of 
square waves 5'^, where i = l,2,3,...,n, do not make up a system of orthogonal 
functions. 
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(5.2) : The dashed line indicates the approximation to /(t), specified in 

i=l 


Figlte 5 




























































































































































































The results obtained upon carrying out the type of analysis described of 
a function characterized in an interval At divided into n sub-intervals with an 
equal duration can be presented in a sequence of dyads (ordered pairs) such that 
the first element of the first dyad is the frequency /i corresponding to and 
the second element of the first dyad is the coefficient Ci ; the first element of the 
second dyad is the frequency /2 corresponding to S 2 and the second element 
of that dyad is the coefficient C 2 ; and so on successively, such that the first 
element of the dyad is the frequency /„ corresponding to Sn and the second 
element of that n*’’ dyad is the coefficient (7„. 

Consider, for example, the sequence of 18 dyads obtained when carrying out 
the type of analysis described of the f{t) specified in Q, if the interval At is 
divided into 18 sub-intervals; 


(/ijCi) = (0.1250000; 117.12980) 
(MCs) = (0.1406250;-210.98830) 
(/5;C'5) = (0.1607143; 9.35088) 
{fr,C7) = (0.1875000; 61.27212) 
(/gjCg) = (0.2250000; 12.81335) 
(/ibC'ii) = (0.2812500;-85.68506) 
(/i3;Ci 3) = (0.3750000; 8.51973) 
(AsiCis) = (0.5625000;-69.86421) 
ifi7;Ci7) = (1.1250000;-9.26140) 


(A; C 2 ) =, (0.1323529; 50.27631) 
(/4; Ci) = (0.1500000; -53.27896) 
(/gjCe) = (0.1730769; 12.58025) 
(/gjCg) = (0.2045455; 49.80105) 
(/io;C'io) = (0.2500000; 4.03101) 
ifi2-,Ci2) = (0.3214286; 12.88482) 
(/ 14 ; Cii) = (0.4500000; 60.38772) 
(/leiC'ie) = (0.7500000; 28.08997) 
ifis-,Cis) = (2.2500000; -32.60758) 


This approximation to function /(t) specified in 0 can be expressed in 
the frequency domain. To achieve this objective, for each of the frequencies 
considered /i, / 2 , /s, • • •, /is, the corresponding coefficients Ci, C 2 , C 3 ,..., C 18 
must be indicated. 

The expression in the frequency domain of this approximation to f{t) will 
be called the Square Wave Transform (SWT) of that approximation to f{t). 
Note that previously (in figuresand 5, for example), each approximation to 
f{t) was represented in the time domain. This SWT is displayed in figure]^ 
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Figure 6: SWT of the approximation to f{t) obtained by dividing At into 18 
sub-intervals. 

Of course, the SWTs corresponding to numbers as large as desired of equal 
sub-intervals into which At is divided can be obtained for the f{t) specified 
in or for any other function of the time which, in a particular interval At, 
satisfies the conditions of Dirichlet. 

In (7.1), (7.2), and (7.3) of figure 7, the SWTs obtained for the approxima¬ 
tions to the /(t) specified in 0 are shown for n = 100, n = 1000, and n = 2000, 
respectively. 
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(7.1) : SWT of the approximation to the /(t), specified in equation (^, for n = 100. 



(7.2) : SWT of the approximation to the /(t), specified in equation 0 , for n = 1000. 

Figure 7 
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(7.3) : SWT of the approximation to the /(t), specified in equation (Q, for n = 2000. 

Figure 7 


In (8.1), (8.2), (8.3), and (8.4) of figure 8, partial representations can be 
seen (up to frequency fi = 2) of the SWTs of the different approximations to 
the f{t) specified in 
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(8.1) : Partial representation of the SWT of the approximation to the f(t) specified in 0 for 
n = 1000 



(8.2) : Partial representation of the SWT of the approximation to the f(t) specified in 0 for 
n = 2000 


Figure 8 
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(8.3) : Partial representation of the SWT of the approximation to the f(t) specified in 0 for 
n = 4000 



(8.4) : Partial representation of the SWT of the approximation to the f(t) specified in 0 for 
n = 8000 


Figure 8 
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In figure 8, it can be observed that in all 4 cases considered the “promi¬ 
nent coefficients” correspond to certain frequencies. (They are prominent in the 
sense that their moduli are quite larger than the moduli of the coefficients cor¬ 
responding to frequencies near to those considered.) They have been indicated 
by the letters A, B, C, and D. The dyads corresponding to the coefficients for 
these cases are the following: 

n = 1000 

A: (0.2441410; 541.50054) 

B: (0.4882812; 270.06817) 

C: (0.9765633; 134.80741) 

D: (1.9531250; 66.39768) 
n = 2000 

A: (0.2441410; 342.97162) 

B: (0.4882812; 171.42003) 

C : (0.9765633; 85.45887) 

D: (1.9531250; 42.22607) 
n = 4000 

A: (0.2441410; 209.36325) 

B: (0.4882812; 104.64873) 

C : (0.9765633; 52.19877) 

D: (1.9531250; 25.84766) 
n = 8000 

A: (0.2441410; 541.40777) 

B: (0.4882812; 62.68744) 

C : (0.9765633; 31.28092) 

D: (1.9531250; 15.51459) 

Note that in these 4 cases, the frequency corresponding to coefficient A is 
equal to 0.2441410; the frequency corresponding to the coefficient B is equal to 
0.4882812; the frequency corresponding to the coefficient C is equal to 0.9765633; 
and the frequency corresponding to the coefficient D is equal to 1.9531250. In 
other words, the frequencies of those “prominent coefficients” are invariable in 
the changes specified for the value of n from 1000 to 2000, from 2000 to 4000, 
and from 4000 to 8000. 
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3 Analysis of Sequences of Samples (Measured 
Values) from an Electromyographic Record¬ 
ing 

The SWM and the corresponding expression of the results obtained with the 
SWT can be used for the analysis of sequences of “samples”, or measured values, 
from different types of recordings. In particular, they can be used for the anal¬ 
ysis of recordings which are important in medicine, as are those of the electro¬ 
cardiogram (ECG), the electroencephalogram (EEG), and the electromyogram 
(EMG). Results obtained by using the SWM for the analysis of a sequence of 
samples from an electroencephalographic recording were described in [2]. 

Suppose that, to take samples from a recording, a sampling frequency fs of 
250 Hz is used. In this case let it be admitted that every one-second lapse (1 s) 
has been divided into 250 sub-intervals of equal duration and that the values of 
the samples “correspond” to the values (mentioned in section 2) of the function 
analyzed, at the midpoints of those different sub-intervals. In general, for any 
numeric value of the fs used, the number of sub-intervals into which the unit of 
time 1 s has been divided is equal to that numeric value. 

Thus if with a specific /s, such as one of 250 Hz, a sequence of samples is 
taken for a certain At, equal to 5 s, for example, to compute the total number 
n of samples taken during that At, the fs must be multiplied by that At: 

n = fs' At (4) 

Eor the values of fs and At specified in the above paragraph, the following 
equality is found: 

1250 = 250s“^ •5s 

The SWM may be used to analyze a function /(t) characterized analytically 
in an interval A = 5 s. If n (the number of sub-intervals into which At is divided) 
is equal to 1250, for example, the sequence of values Ei, ^ 2 , V 3 ,..., E 1250 will be 
obtained for the midpoints of the specific sub-intervals as follows: Vi will be 
the value of /(t) at the midpoint of the first sub-interval, V 2 will be the value 
of /(t) at the midpoint of the second sub-interval, and so on successively, with 
Ei 25 o fhe value of /(t) for the last of the sub-intervals considered. 

When operating with fs = 250 Hz, and taking samples for At = 5 s, 1250 
samples will be obtained. To apply the SWM to this sequence of samples, it 
is considered that the first “corresponds” to Ei, the second to E 2 , and so on, 
successively such that the last sample “corresponds” to Ei 25 o- In other words, 
the sequence of samples is treated the same as the sequence Ei, E 2 , E 3 ,..., E 1250 , 
when applying the SWM to an analytically characterized function. 

Equation which makes it possible to compute each fi corresponding to 
each Si (where i = 1,2,3,...,n), remains valid for sequences of samples of a 
recording: 


27 



( 5 ) 


If Q is substituted in ([^, the following equation]^ is obtained: 


2 (/,At-{z-l))’ 


i = 1, 2, 3,..., n 


Note the equations obtained if in a) i = 1; and h) i = n. 

For a), (i = 1) —)► fi = In this first case, fi does not depend on /g, 
but rather only on At. 

For b), (i = n) —^ fn = In this second case, fn does not depend on At, 
but rather only on fs. 

The graph corresponding to an electromyographic recording made during 
At = 5 s with fs = 250 Hz is presented in figure [8] . 



Figure 9: Electromyographic recording made during At = 5s with fs = 250Hz. 

Extending only to /^ = 2 Hz, the graph in figure is a section of the SWT 
corresponding to the recording represented in figure 
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Figure 10; Partial plot of the SWT corresponding to the electromyographic 
recording in figure]^ The moduli of the coefficients not represented (the Ci for 
which fi > 2) are much smaller than the moduli of Ci which were represented. 


The sequence of dyads for values of /j < 0.12 is as follows: 


(fi,Ci) = (0.1000000,1.655) 
(/3,C'3) = (0.1001603,-0.152) 
(/5,C'5) = (0.1003210,0.025) 
(/7,C7) = (0.1004823,0.017) 
(/9,C9) = (0.1006441,0.019) 
(fii,Cn) = (0.1008065,-0.033) 
(fi3,Ci3) = (0.1009693,0.007) 
(/i5,Ci5) = (0.1011327,0.023) 
(fi7,Ci7) = (0.1012966,0.004) 
(/i9,Ci9) = (0.1014610,-0.067) 
(/21,C2i) = (0.1016260,-0.010) 
(f 23 ,C 23 ) = (0.1017915,0.016) 
(f 25 ,C 25 ) = (0.1019576,0.019) 
(/27,C27) = (0.1021242,-0.009) 
(f 29 ,C 29 ) = (0.1022913,-0.021) 


(f2,C2) = (0.1000801,0.002) 
(/4,C4) = (0.1002406,-0.013) 
(f6,Ce) = (0.1004016,0.004) 
(fs,C8) = (0.1005632,-0.026) 
(fio,Cio) = (0.1007252,-0.027) 
(/i2,C'i2) = (0.1008878,-0.021) 
(/i4,Ci4) = (0.1010509,0.011) 
(/i6,Ci6) = (0.1012146,0.012) 
(fi8,Ci8) = (0.1013788,0.006) 
(/20,C2o) = (0.1015435,0.011) 
(/22,C22) = (0.10 17087,0.011) 
(/ 24 ,C 24 ) = (0.1018745,-0.007) 
(/26,C26) = (0.1020408,-0.005) 
(/28,C28) = (0.1022077,0.005) 
(/30,C3o) = (0.1023751,0.028) 
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(/31,C3i) = (0.1024590,0.010) 
(/ 33 ,C 33 ) = (0.1026273,-0.001) 
(/ 35 ,C 35 ) = (0.1027961,0.340) 
(/ 37 ,C 37 ) = (0.1029654,0.006) 
(/ 39 ,C 39 ) = (0.1031353,-0.001) 
(/ 4 i,C 4 i) = (0.1033058,0.017) 
(/ 43 ,C 43 ) = (0.1034768,0.006) 
(/ 45 ,C 45 ) = (0.1036484,0.010) 
(/ 47 ,C 47 ) = (0.1038206,-0.031) 
(/49,C'49) = (0.1039933,0.006) 
(/ 5 i,C 5 i) = (0.1041667, -0.020) 
(/ 53 ,C 53 ) = (0.1043406,0.005) 
(/ 55 ,C 55 ) = (0.1045151,-0.010) 
(/ 57 ,C 57 ) = (0.1046901,-0.012) 
(759,^59) = (0.1048658,0.027) 
ifeuCei) = (0.1050420,0.014) 
{h3,Ces) = (0.1052189,0.020) 
(/65,C65) = (0.1053963,0.006) 
(/67,C67) = (0.1055743,0.018) 
(/69,C69) = (0.1057530,0.003) 
(/ 7 i,C' 7 i) = (0.1059322,-0.034) 
(773,^73) = (0.1061121,-0.025) 
(/ 75 ,C 75 ) = (0.1062925,0.055) 
{fn,Cn) = (0.1064736,0.004) 
(/ 79 ,C 79 ) = (0.1066553,0.018) 
(/81,C8i) = (0.1068376,-0.004) 
(783,^83) = (0.1070205,0.021) 
(785,^85) = (0.1072041,-0.038) 
(787,^87) = (0.1073883,0.002) 
(789,^89) = (0.1075731,0.005) 
(791, C 91 ) = (0.1077586,-0.012) 
( 793 ,^ 93 ) = (0.1079447,0.002) 
(795,^95) = (0.1081315,-0.058) 
(797,^97) = (0.1083189,0.000) 
( 799 ,^ 99 ) = (0.1085069,-0.300) 
(7ioi,Cioi) = (0.1086957,0.014) 
(7 io3,C'io3) = (0.1088850,-0.037) 


(732,^32) = (0.1025431,-0.005) 
( 734 ,^ 34 ) = (0.1027116,0.010) 
(736,^36) = (0.1028807,0.024) 
(738,^38) = (0.1030503,0.016) 
(740, C40) = (0.1032205,0.001) 
(742, C42) = (0.1033912,0.006) 
( 744 , 6 * 44 ) = (0.1035626,-0.019) 
(746,^46) = (0.1037344,-0.015) 
(748,6*48) = (0.1039069,-0.011) 
(750,^50) = (0.1040799,0.015) 
(752,6^52) = (0.1042535,-0.031) 
(754,^54) = (0.1044277,0.011) 
(756,^56) = (0.1046025,-0.004) 
(758,^58) = (0.1047779,0.010) 
(760,6*60) = (0.1049538,-0.005) 
(762,^62) = (0.1051304,-0.026) 
(764,6*64) = (0.1053075,0.013) 
(766,^66) = (0.1054852,0.010) 
(768,^68) = (0.1056636,-0.002) 
(770, C70) = (0.1058425,0.011) 
( 772 , 6 * 72 ) = (0.1060221,-0.005) 
( 774 , 6 * 74 ) = (0.1062022,-0.039) 
(776,6*76) = (0.1063830,-0.005) 
(778,^78) = (0.1065644,0.034) 
(780,^80) = (0.1067464,0.010) 
(782,^82) = (0.1069290,-0.011) 
(784,^84) = (0.1071123,-0.010) 
(786,^86) = (0.1072961,-0.009) 
(788,^88) = (0.1074807,0.013) 
(790,6*90) = (0.1076658,-0.009) 
( 792 , 6 ^ 92 ) = (0.1078516,0.007) 
( 794 , 6 * 94 ) = (0.1080380,0.003) 
(796,^96) = (0.1082251,-0.010) 
(798,^98) = (0.1084128,0.013) 
(7ioo,Cioo) = (0.1086012,0.013) 
(7 io2,Cio2) = (0.1087903,-0.013) 
(7104,6*104) = (0.1089799,-0.008) 
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(/io5,Cio5) = (0.1090750,0.001) 
(/io7,Cio7) = (0.1092657, -0.146) 
(/io9,C'io 9) = (0.1094571,0.019) 
(/iii,C'ni) = (0.1096491,-0.014) 
(/ii3,Cii3) = (0.1098418,0.005) 
(/ii5,Cii5) = (0.1100352,0.021) 
(/ii7,C'ii 7) = (0.1102293,-0.037) 
(/ii 9 ,C'n 9 ) = (0.1104240,0.056) 
(/i21,Ci2i) = (0.1106195,0.003) 
(/i23,Ci23) = (0.1108156,0.007) 
(/i25,C'i25) = (0.1110124,-0.009) 
(/i27,C'i27) = (0.1112100,0.019) 
(/i29,Ci29) = (0.1114082,0.043) 
(/i3i,Ci3i) = (0.1116071,0.128) 
{fi3Z,Cn3) = (0.1118068,-0.042) 
(/i35,Ci35) = (0.1120072,0.000) 
(/i37,Ci37) = (0.1122083,0.015) 
(/i39,C'i39) = (0.1124101,-0.010) 
(/i41,C'i4i) = (0.1126126,-0.005) 
(/i43,Ci43) = (0.1128159,0.035) 
(/i45,Ci 45) = (0.1130199,0.007) 
(/i 47 , C' 147 ) = (0.1132246, -0.024) 
(/i49,C'i49) = (0.1134301,-0.020) 
(/i51,Ci5i) = (0.1136364,-0.008) 
(/i53,Ci53) = (0.1138434,0.015) 
(/i55,C'i55) = (0.1140511,-0.014) 
(/i57,C'i57) = (0.1142596,0.001) 
(/i59,Ci59) = (0.1144689,-0.002) 
(/i6i,Ci6i) = (0.1146789,-0.014) 
(/i63,C'i63) = (0.1148897,0.824) 
(/i65,Ci65) = (0.1151013,-0.003) 
(/i67,Ci67) = (0.1153137,0.010) 
(/i69,C'i69) = (0.1155268,0.005) 
(/i71,C'i7i) = (0.1157407, -0.044) 
(/i73,Ci73) = (0.1159555,-0.012) 
(/i75,Ci75) = (0.1161710,0.005) 
(/i77,C'i77) = (0.1163873,0.004) 


{fw6,Cwe) = (0.1091703,-0.005) 
{fios,Cws) = (0.1093613,0.000) 
{fiio,Ciio) = (0.1095530,-0.022) 
(/ii 2 ,C'n 2 ) = (0.1097454,-0.034) 
(/ii 4 ,Cii 4 ) = (0.1099384,0.003) 
(/ii 6 ,Cii 6 ) = (0.1101322,0.019) 

= (0.1103266,-0.023) 
(/i 20 ,C'i 2 o) = (0.1105217,0.019) 
(/i 22 ,Ci 22 ) = (0.1107174,0.008) 
(/i24, Cm) = (0.1109139,-0.009) 
(/i26,Ci26) = (0.1111111,0.002) 
(/i28,Ci28) = (0.1113090,0.016) 
(/i30,Ci3o) = (0.1115076,-0.008) 
(/i32,Ci32) = (0.1117069,-0.015) 
(/i34,Ci34) = (0.1119069,0.019) 
(/i36,Ci36) = (0.1121076,-0.008) 
(/i38,Ci38) = (0.1123091,0.008) 
(/i40,Ci4o) = (0.1125113,-0.023) 
(/i42,Ci42) = (0.1127142,0.009) 
(/i44,Ci44) = (0.1129178,-0.002) 
(/i46,Ci46) = (0.1131222,0.032) 
(/i48,Ci48) = (0.1133273,-0.002) 
(/i50,Ci5o) = (0.1135332,0.000) 
(/i52,Ci52) = (0.1137398,0.006) 
(/i54,Ci54) = (0.1139471,0.011) 
(/i56,Ci56) = (0.1141553,0.014) 
(/i58,Ci58) = (0.1143641,0.021) 
(/i60,Ci6o) = (0.1145738,0.013) 
(/i62,Ci62) = (0.1147842,0.005) 
(/i64,Ci64) = (0.1149954,0.012) 
(/i66,Ci66) = (0.1152074,0.005) 
(/i68,Ci68) = (0.1154201,0.005) 
(/i70,Ci7o) = (0.1156337,-0.005) 
(/i72,Ci72) = (0.1158480,-0.021) 
(/i 74 ,Ci 74 ) = (0.1160631,-0.008) 
(/i76,Ci76) = (0.1162791,-0.001) 
(/i78,Ci78) = (0.1164958,-0.018) 


31 



(/i79,Ci79) = (0.1166045,0.097) 
(/i8i,Ci8i) = (0.1168224,-0.004) 
(/i 83 ,C'i 83 ) = (0.1170412,0.017) 
(/i 85 ,C'i 85 ) = (0.1172608,-0.017) 
(/i 87 ,Ci 87 ) = (0.1174812,-0.013) 
(/i 89 ,Ci 89 ) = (0.1177024,0.009) 

= (0.1179245,-0.027) 
(/i93,C'i93) = (0.1181474,0.034) 
(/i95,Ci95) = (0.1183712,-0.404) 
(/i97,Ci97) = (0.1185958,-0.044) 
(/i99,C'i99) = (0.1188213,0.032) 
(/20i,C'2oi) = (0.1190476,0.018) 
(/ 203 ,C 203 ) = (0.1192748,0.025) 
(/ 205 ,C 205 ) = (0.1195029,0.006) 
(/ 207 ,C' 207 ) = (0.1197318,0.024) 
(/ 209 ,C 209 ) = (0.1199616,0.005) 


(/i80,Ci8o) = (0.1167134,-0.004) 
(/i82,Ci82) = (0.1169317,0.014) 
(/i84,C'i84) = (0.1171509,-0.032) 
(/i86,C'i86) = (0.1173709,0.003) 
(/i88,Ci88) = (0.1175917,0.010) 
(/i90,Ci9o) = (0.1178134,0.008) 
(/i92,C'i92) = (0.1180359,-0.009) 
(/i94,C'i94) = (0.1182592,0.000) 
(/i96,Ci96) = (0.1184834,-0.009) 
(/i98,Ci98) = (0.1187085,-0.023) 
(/20o,C'2oo) = (0.1189343,0.007) 
(/202,C'202) = (0.119 1 611,0.011) 
(/ 204 ,C 204 ) = (0.1193887,-0.008) 
(/ 206 ,C 206 ) = (0.1196172,0.032) 
(/ 208 ,C' 208 ) = (0.1198466,-0.012) 


In 0 potential users can find a tool which makes it possible to obtain the 
SWT of sequences of samples of different recordings, for n < 3000. 


4 Analysis of an Image 

The SWM is also useful for the analysis of images and the corresponding results 
can be expressed with the SWM. 

First, a “pseudo-image” (shown in figure pT]) has been chosen to illustrate 
how the SWM may be used for analysis purposes. The method applied for the 
analysis of a pseudo-image is the same as that used for the analysis of genuine 
images. The term “pseudo-image” indicates that no reference is made to a 
genuine image of something real, such as an object or a living being, or abstract 
art with aesthetic value. Given that the method used, the SWM, for the analysis 
of a genuine image and of a pseudo-image is the same, the latter will serve to 
provide a clear and simple example of how that method can be applied. 
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Figure 11: Representation of a pseudo-image. The name of each pixel is given 
at the top, and the level of gray corresponding to that pixel is given at the 
bottom. 

The name of each pixel in the pseudo-image displayed in figure is given 
at the top. The numerical value (ranging from 0 for black to 255 for white) 
corresponding to the level of gray of the pixel is given at the bottom. (The 
levels of gray were chosen arbitrarily.) 

Although figure[^can be considered as a representation of a matrix of pixels, 
the diverse elements in that matrix (Pij, where i = 1, 2, 3,4; and j = 1, 2, 3,4) 
are referred to differently from what is usually found. Here it has been taken 
into account that in the system of orthogonal Cartesian coordinates used (in 
which the x-axis is that of the abscissas, and the ^-axis is that of the ordinates), 
the point (i,j) (where i = 1,2,3,4; and j = 1,2, 3,4) whose abscissa is i and 
whose ordinate is j, is the center point of the pixel Pij. 

In figure it has been shown how each of the axes {x and y) in figure 11 
has been processed with the same approach as that of figure with the x-axis 
(the only axis of coordinates used in this case). 
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Figure 12: How to apply the SWM to the analysis of the pseudo-image repre¬ 
sented in figure 11. (See indications in text.) The superscript 1 located on the 
left of the coefficient indicates that the coefficient corresponds to the x-axis and 
the superscript 2 on the left of the coefficient indicates that it corresponds to 
the ^-axis. 

Note, for example, how the equation corresponding to pixel P 4^3 is formed: 
that pixel belongs to both column 4 and row 3 in figure 

1. Take the sequence of these two sets: a) the set of elements underneath 
column 4 in figure 12; and b) the set of elements at the left of row 3 in 
the same figure: 

('Cl, -1^2, -'C3, -IC4) 2^2, 2C4) 


Note that the above expression is an ordered pair of sets. 

2. The Cartesian product is found for the two sets corresponding to that 
ordered pair of sets: 

(1^ 1/^ 1/^ \ s. (2/^ 2^ 2^ 2r^ \ _ p* _ 

V ^25 ^35 '^4j ^ V ^1’ ^25 ^35 ~ ^4,3 ~ 

{( 'Cl, ^C,), ( 2 ^ 2 ), ( ( 'Cl, 2 ^ 4 ), 

(-'C2,2Ci), (-'C2,2C2), (-'C2,-2C3), (-'C2,2C4), 
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(-V3,^Ci), {-^c„^c,), (-^Cs,-^Cs), i-^Cs,^C,), 

{-^C„^C,), {-^C„^C,), {-^c„-^c,), (-^C„^C,)} 

The preceding Cartesian product corresponding to pixel has been 
designated as P 4 3 . 

3. The set P 4 3 is acted on by an operator O such that: 

(a) it converts each element P 4 3 (a certain ordered pair) into a single 
element which will be a coefficient with two subscripts: the first is 
the same as that of the first element of the ordered pair, and the 
second subscript is the same as that of the second element of the 
ordered pair; and 

(b) if the signs preceding the two elements of the ordered pair have the 
same sign (+ and +; or — and —), the coefficient is positive; if they are 
different (+ and —; or — and +), the coefficient obtained is negative. 

Therefore, for the case considered, the following result is obtained: 

0(P4*3)=(C'i,1,C^1,2,-Ci,3,Ci,4, 

— ^ 2 , 1 ? —^ 2 , 2 ? ^2,3? —^2,4? 

— Cs,!^ ~C'3,2, ^3,3, —C3,4, 

— C'4p, —C 4 ^ 2 , C'4^3, —04^4) 

The right-hand member of the above equation is a set of coefficients, each 
of which has two subscripts and is preceded by the + sign or by the — 
sign. That set will be called Of course, ^ must not be confused 

with the coefficient (^ 4 ^ 3 . 

4. The algebraic sum is figured for all the elements in and the result is 
equated to the numeric value of the level of gray corresponding to pixel 
P 4,3 (he., 255). 

Therefore, the linear algebraic equation corresponding to pixel P 4^3 is ob¬ 
tained: 


C^l,! + ^71^2 — h7i^3 + Ci^4 — C2,l — C2,2 + C2,3 ~ ^^2^4 

— Cs,! — Cs^2 + ^^ 3^3 — C3,4 — C4q, —C4,2 + ^^ 4^3 — C4,4 = 255 

If the same procedure is applied to the set of elements under column 3 in figure 
12 and to the set of elements appearing at the left of row 1 in that same figure, 
the following linear algebraic equation may be specified for pixel P 3 q: 

^1,1 ff Ci^2 ff h7i^3 + Ci^4 + h72,l H“ ^^ 2^2 ff ^2^3 H“ ^^ 2^4 

~ — C3^2 ~ — (^3^4 -h C4^i “h ^ 4^2 + ^^ 4^3 (^ 4^4 = 25 
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If the same procedure is carried out for each of the remaining pixels of the 
pseudo-image in figure 11, the following system of linear algebraic equations is 
obtained: 


^1,1 

+ C'l,2 

+ Ci,s 

+ Ci,4 

+ C2,l 

+ C2,2 

+ C2,S 

+ C2,4 




+ C's,! 

+ Cs,2 

+ Cs,s 

+ Cs,4 

+ C4,i 

+ C4,2 

+ C4,s 

+ C4,4 

= 100 

Cl.l 

+ C'l,2 

+ Ci,s 

-Ci,4 

+ C2,l 

+ C2,2 

+ C2,S 

-C2,4 




+ C'sq 

+ Cs ,2 

+ Cs,s 

— Cs,4 

+ C4,i 

+ C4,2 

+ C4,s 

-C4,4 

= 98 


+ C'l,2 

-Ci,s 

+ Ci,4 

+ C2,l 

+ C2,2 

— C2,S 

+ C2,4 




+ C's,! 

+ Cs,2 

-Cs,s 

+ Cs,4 

+ C4,i 

+ C4,2 

-C4,S 

+ C4,4 

= 195 

C'l,! 

-Cl,2 

-Ci,s 

-Ci,4 

+ C2,l 

-C2,2 

— C2,S 

-C2,4 




+ Cs,! 

— Cs,2 

— Cs,s 

— Cs,4 

+ C4,i 

— C4,2 

— C4,s 

-C4,4 

= 55 


+ Cl,2 

+ Ci,s 

+ Ci,4 

+ C2,l 

+ C2,2 

+ C2,S 

+ C2,4 




+ Cs,! 

+ Cs,2 

+ Cs,s 

+ Cs,4 

-C4,l 

-C4,2 

— C4,s 

-C4,4 

= 38 


+ Cl,2 

+ Ci,s 

-Ci,4 

+ C2,l 

+ C2,2 

+ C2,S 

-C2,4 




+ Cs,! 

+ 

to 

+ Cs,s 

-Cs,4 

-C4,l 

-C4,2 

-C4,S 

+ C4,4 

= 3 


+ Cl,2 

-Cl,3 

+ Ci,4 

+ C2,l 

+ C2,2 

— C2,S 

+ C2,4 




+ Cs,! 

+ Cs ,2 

-Cs,s 

+ Cs,4 

-C4,l 

— C4,2 

+ C4,s 

-C4,4 

= 6 

Ci,i 

-Ci,2 

— Ci,s 

-Ci,4 

+ C2,l 

— C2,2 

— C2,S 

-C2,4 




+ Cs,i 

— Cs ,2 

— Cs,s 

-Cs,4 

— C4,i 

+ C4,2 

+ C4,s 

+ C4,4 

= 4 

ClA 

+ Cl,2 

+ Ci,s 

+ Ci,4 

+ C2,l 

+ C2,2 

+ C2,S 

+ C2,4 




-Cs,i 

— Cs ,2 

-Cs,s 

— Cs,4 

+ C4,i 

+ C4,2 

+ C4,s 

+ C4,4 

= 25 

Ci,i 

+ Cl,2 

+ Ci,s 

-Ci,4 

+ C2,l 

+ C2,2 

+ C2,S 

-C2,4 




-Cs,i 

— Cs ,2 

— Cs,s 

+ Cs,4 

+ C4,i 

+ C4,2 

+ C4,S 

— C4,4 

= 77 

ClA 

+ Cl,2 

-Ci,s 

+ Ci,4 

+ C2,l 

+ C2,2 

-C2,S 

+ C2,4 




-Cs,i 

— Cs ,2 

+ Cs,s 

— Cs,4 

+ C4,i 

+ C4,2 

-C4,S 

+ C4,4 

= 249 

c*!,! 

-Ci,2 

— Ci,s 

-Ci,4 

+ C2,l 

— C2,2 

— C2,S 

-C2,4 




-Cs,i 

+ Cs ,2 

+ Cs,s 

+ Cs,4 

+ C4,i 

-C4,2 

— C4,s 

-C4,4 

= 69 


+ Cl,2 

+ Ci,s 

+ Ci,4 

-C2,l 

— C2,2 

-C2,S 

— C2,4 




-Cs,i 

— Cs ,2 

-Cs,s 

— Cs,4 

-C4,l 

— C4,2 

-C4,S 

-C4,4 

= 55 


+ Cl,2 

+ Ci,s 

-Ci,4 

-C2,l 

— C2,2 

— C2,S 

+ C2,4 




-Cs,i 

— Cs ,2 

— Cs,s 

+ Cs,4 

— C4,i 

-C4,2 

-C4,S 

+ C4,4 

= 12 

C*!,! 

+ Cl,2 

-Ci,s 

+ Ci,4 

— C2,l 

-C2,2 

+ C2,S 

-C2,4 




-Cs,i 

— Cs ,2 

+ Cs,s 

— Cs,4 

-C4,l 

— C4,2 

+ C4,s 

-C4,4 

= 255 


-Ci,2 

-Ci,s 

— Ci,4 

-C2,l 

+ C2,2 

+ C2,S 

+ C2,4 




-Cs,i 

+ Cs ,2 

+ Cs,s 

+ Cs,4 

-C4,l 

+ C4,2 

+ C4,s 

+ C4,4 

= 81 


The 16 unknowns of the preceding system of linear algebraic equations are 
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the 16 coefficients where i = 1, 2, 3,4; and j = 1, 2, 3,4. 

If the system of equations is solved, the following values are found for 
the 16 unknowns: 


Ci,i = 112.50000 
Cl,2 = -78.50000 
Cl,3 = 15.25000 
Cl,4 = 28.25000 
C2,i = 27.50000 
C2,2 = 22.25000 
C2,3 = -23.50000 
C2,4 = 42.75000 


C3,i = -34.00000 
C3,2 = -14.00000 
C3,3 = 32.25000 
C3,4 = -31.75000 
C4,i = 51.00000 
C4,2 = -55.25000 
C4,3 = 13.50000 
C4,4 = -8.25000 


Information will be provided later about the analysis of the image of Lenna 
mentioned above, regarding the meaning of this type of coefficients. For now, it 
will suffice to say that the system of equations (6) obtained for the pseudo-image 
considered “corresponds” to the system of equations (3) obtained by analyzing 
the function of the single variable specified in (1). 

In general, to determine the equation corresponding to any pixel Pij (where 
i = 1, 2, 3,..., n; and j = 1, 2, 3,..., n) of an image composed of n rows and n 
columns of pixels (i.e., with pixels), the following procedure may be used: 


1. Using the same criterion as that of the case of the pseudo-image discussed 
above, it is possible to determine the set of elements underneath column i, 
and the set of elements at the left of row j of the image analyzed. Recall 
that the abscissa i of the center of pixel Pij indicates the column where 
that pixel is located; and j, the ordinate of the center of Pij indicates the 
row in which that pixel is found. 

2. The Cartesian product is found for these two sets, and the result obtained 
is referred to as P*^-. 

3. An action is carried out on the set P*^ by an operator O such that: 

(a) it converts each element of P*^- (a particular ordered pair) into a 
unique element which will be a coefficient with two subscripts: the 
first is the same as the subscript of the first element of that ordered 
pair and the second subscript is the same as the subscript of the 
second element of that ordered pair; and 

(b) if the signs preceding the two elements of the ordered pair are the 
same (both positive or both negative), then the sign of the coefficient 
obtained is positive. If, on the other hand, the signs are different 
(positive and negative, or negative and positive), then the sign of the 
coefficients is negative. 
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4. C^j is used to refer to the set obtained by having the operator O act on 
the set P*-: 



5. The algebraic sum of all the elements is of is calculated and the result 
is equated to the numeric value of the level of gray corresponding to pixel 
P 

When doing the same with each pixel P^j (where i = 1,2,3, and 

j = 1, 2, 3,..., n), the result is a system of m? linear algebraic equations with 
unknowns. Those unknowns, which can be found, are the coefficients Ci^j 
(where i = 1, 2, 3,..., n; and j = 1, 2, 3,..., n). 

The classic image of Lenna to be analyzed using the SWM is displayed in 
figure 13. 


y 



X 


Figure 13: Image of Lenna 


The image in figure is composed of 262,144 pixels. It can be considered 
a matrix of pixels with 512 rows, each made up of 512 pixels, and 512 columns, 
each of which is also made up of 512 pixels. Each of the 262,144 pixels belongs 
to only one of the 512 rows and only one of the 512 columns. In other words, 
the image of Lenna analyzed can be viewed as a matrix arranged in a 512 x 512 
square grid of pixels. 
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If one applies the SWM directly to that image of Lenna, a system of 262,144 
linear algebraic equations is obtained. Because the authors of this paper cannot 
solve this system of linear equations due to the limitations of the computational 
tools available, the image has been further divided into 256 sub-images, each of 
which is a 32 X 32 matrix, consisting of: 

1. 16 rows of sub-images, each composed of 16 sub-images; or 

2. 16 columns of sub-images, each composed of 16 sub-images. 

In other words, the image of Lenna analyzed can also be considered as a 
matrix arranged in a 16 x 16 square grid of sub-images. 

The computational tools available to the authors do have the capacity to 
apply the SWM to the analysis of each of these 256 sub-images, and this process 
was carried out successfully. 

The columns of sub-images have been numbered from left to right, from the 
first column (1) to the last column (16). The rows of sub-images have also 
been numbered from bottom up, row 1 to row 16. The sub-image belonging to 
column k (/c = 1, 2, 3,..., 16) and row I (/ = 1, 2, 3,..., 16) is called 

In figure 14, it can be seen how the image of Lenna has been divided into 
256 sub-images, so that each sub-image belongs only to one column and one 
row of the matrix. 
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The notations introduced above for the pixels and for the sub-images can 
be used to refer to any pixel of the image of Lenna, by using the letter “P” for 
pixel and a tetrad of subscripts of that letter. Therefore, P 4 , 17 , 8,14 refers to the 
pixel located in the fourth column and the seventeenth row of the sub-image 
located, in turn, in the eighth column and the fourteenth row of the matrix of 
sub-images of Lenna. 

The analysis of sub-image of the image of Lenna (i.e., the sub-image 
located in column 9, row 8 of the matrix of the sub-images of the image analyzed) 
is discussed below. The same approach was used for the analysis of all 255 
remaining sub-images of the image of Lenna. Sub-image Jg^g of the image of 
Lenna is displayed in figure 15. 

y 


X 

Figure 15: Sub-image /g,8 of the image of Lenna 

When using the SWM to analyze /g^8 5 a system of 1024 linear algebraic 
equations was solved. Recall that each sub-image was composed of a matrix of 
32 X 32 pixels. For each of those 1024 pixels, a linear algebraic equation was 
obtained using the procedure described. 

The left-hand member of the linear algebraic equation obtained for each of 
the 1024 pixels of /g^8 is an algebraic sum of the 1024 coefficients Ci^j (where 
i = 1, 2, 3,..., 32; and j = 1, 2, 3,..., 32); and the right-hand member of the 
equation is the numeric value of the level of gray corresponding to the pixel 
considered. When solving this system of 1024 equations, each of the 1024 coef¬ 
ficients (the unknowns of the system of equations) can be positive, negative or 
null. Hence, according to the algebraic equation corresponding to each pixel, 
each coefficient can make a positive, negative or null contribution to that alge¬ 
braic sum. Of course, to determine the contribution of each coefficient to that 
algebraic sum, both the computed value of each coefficient and the preceding 
sign in the algebraic sum must be taken into account. Thus, for example, when 
solving the system of equations, if it is determined that the value of a coefficient 
is —250, and the sign preceding it is negative, the contribution of that coeffi¬ 
cient to that algebraic sum is equal to a positive numeric value: 250. However, 
if the sign preceding that coefficient is positive, then the contribution of that 
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coefficient to the algebraic sum is negative: —250. 

Figure 16 represents the type of contribution (either positive or negative) 
which the coefficient Ci 7,25 = 14.500 gives to the value of gray of each pixel 
of Ig^s^ according to its numeric value (which then remains unchanged for all 
of the pixels of Ig^g) and to the preceding sign, which can change according to 
the algebraic equation corresponding to each pixel of Ig^g. Blue is used for the 
positive contributions of Ci 7,25 5 and red for the negative contributions for each 
pixel of Ig^s- la figure 16, consideration was not given to the absolute value of 
these contributions to the numeric value of the level of gray of each pixel; it is 
only shown whether they are positive or negative. 



Figure 16: Type of contribution (positive or negative) of Ci 7,25 to each pixel 
Ig^s- The positive contributions are represented in blue; the negative, in red. 


Like figure 16, figure 17 indicates whether the coefficient (^ 29,1 provides a 
positive or negative contribution to each pixel Jg^g- 
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Figure 17: Type of contribution (positive or negative) of C 29 ,i to each pixel of 
The positive contributions are represented in blue; and the negative, in 

red. 

The contributions of each coefficient Cij (where i = 1, 2, 3,, 32; and j = 

1, 2, 3,..., 32) to each pixel of have certain bidimensional patterns. If the 
modulus of each contribution is taken into account in addition to its positive or 
negative sign, these bidimensional patterns are the elements “corresponding” to 
the trains of square waves Si^ S 2 , • - •, Sn- 

An image of the type analyzed can be considered as a function of two vari¬ 
ables: X and y. In this case, to obtain the different approximations to the image 
analyzed, one must add, for each pixel in every sub-image, the patterns men¬ 
tioned, taking into account the modulus of the coefficient corresponding to each 
pattern (just as in the cases of the functions of one variable, the trains of square 
waves are added to obtain the approximations to the functions analyzed). 

In the case of a one-variable function, each coefficient Ci (where i = 1,2,3,...,n) 
is associated to a specific frequency fi. Given the way in which the coefficients 
Cij (where i = 1, 2, 3,..., n; and j = 1, 2, 3,..., n) have been characterized, it is 
clear that each of the latter coefficients is similarly associated to the two spatial 
frequencies fi^x and corresponding to the x-axis and ^-axis respectively. 
The first of these frequencies fi^x characterizes the periodicity made evident by 
the pattern corresponding to each coefficient Cij, according to the x-axis, and 
the second of these frequencies fj^y characterizes the periodicity shown by the 
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same pattern, according to the ^-axis. 

The computations of fi^x and fj^y are carried out by using the following 
equations (7) and (8), obtained from equation (2), substituting At by Ax and 


Ay, respectively (given that 
n = 32: 


and fj^y are spatial frequencies) and equating 


fi,x = 




1 


2Ax 

1 


32 


32 


- (^■ 
32 


1 ) 


i = 1,2,3,..., 32 


j — 1, 2,3,..., 32 


(7) 


( 8 ) 


2Ay \32-{j-l) 

In this case, rather than taking the length of the side of a pixel as the unit 
of length. Ax (which has the same length as Ay) can be used for that purpose. 
Both the length of Ax and that of Ay are equal to the product of 32 times 
the length of one side of a pixel. If that is done, the following equations are 
obtained for fi^x and fj^y'. 


1 

( V 

i = 1,2,3,..., 32 

(9) 

2 ' 

132-(z-1);’ 

ii 

^ 32 \ 

j = 1, 2,3,..., 32 

(10) 

2 ' 

v32-(i-l)j’ 


In the first place, equations and ( P!q| ) are used to consider the case pre¬ 
sented in figure 16, which represents the contribution of (7i7^25 (either positive 


or negative) to each pixel of Jg^g- Given that in this case i = 17 and j = 25, the 
following equations are obtained: 


fi,x — ^ 


32 


2 V32-(17-l)y 2 V32-16) 


1 


32 


2 vie; 


fj,y 


32 


32 - (25 - 1) 


1 / 32 \ _ 1 p2\ _ 4 _ ^ 

2 V32 - 24) J ~ 2\Y) ~ 2~ 


The meaning of the first of the preceding equations is as follows: In interval 
Ax there is one wave on the x-axis, of the pattern represented in figure(See 
figure to confirm.) 

The meaning of the second of the preceding equations is as follows: In inter¬ 
val Ay there are two waves on the ^-axis, of the pattern represented in figure 
(See figure to confirm.) 

Secondly, equations ^ and (10) will be used to consider the case in figure pT| 
which represents the contribution of C 29 ,i (indicating either positive or negative) 
to each pixel Jg^g- Given that in this case i = 29 and j = 1, the following 
equations are obtained: 


fi,X — rj 


32 


2 V32 -(29 - 1) 


32 


32 - 28) 


2 V 4 ; 2 
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fjiV ~ o 


32 


2 V32-(l-l)y 2 V32-0) 


1 


32 


1 _ 1 
2 2 


The meaning of the first of the preceding equations is as follows: In interval 
Ax, there are 4 waves on the x-axis, of the pattern represented in figure 17. 
(See figure to confirm.) 

The meaning of the second of the preceding equations is as follows: In the 
interval A^, there is a half-wave on the ^-axis, of the pattern represented in 
figure 17. (See figure to confirm.) 

In the case of the images, the SWM generates two spatial frequencies, fi^x 
and fj^y, associated with Cij, where i = 1, 2, 3,..., n, and j = 1, 2, 3,..., n. 
Hence, when applying the method it is also possible in the frequency domain to 
represent the results obtained. This expression in the frequency domain of the 
results obtained when analyzing images with SWM will be known as SWT, just 
as when functions of one variable were analyzed. The SWT will be used here 
to express some of the results produced. 

If one adds, for each pixel of every sub-image of Lenna, the positive (or 
negative) contributions corresponding to the different coefficients, one obtains 
once again the image analyzed of Lenna. The analysis process carried out with 
the SWM, as indicated above, has made it possible to determine, unambiguously, 
for each pixel of the sub-image of Lenna, which patterns must be added in 
order to obtain the image of Lenna again. (Those patterns, as shown above, 
“correspond” to the trains of square waves obtained when analyzing functions 
with one variable with the SWM.) 

Suppose that not all but only part of the patterns considered are added. 
For instance, let us admit that for each pixel of every sub-image of Lenna, the 
contributions corresponding to the coefficients Cij such that i < S and j < S 
are the only ones added. The image thus obtained will be called “approximation 
8_32” to the image of Lenna. On the other hand, if for each pixel of every sub¬ 
image of Lenna, the contributions corresponding to the coefficients Cij such 
that i < 2b and j < 25 are the only ones added, the image thus obtained will 
be referred to as “approximation 25_32” to the image of Lenna. Of course, if for 
each pixel of every sub-image of Lenna, the contributions corresponding to all 
the coefficients Cij are added, the “approximation 32_32”, which is the same as 
the image itself, will be produced. 

Several approximations of the image of Lenna have been displayed in |(18.1)| 
to |(18.5)| 


44 
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(18.1) : Approximation 8_32 to image of Lenna 



(18.2) : Approximation 16_32 to image of Lenna 

Figure 18 
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(18.3) : Approximation 24_32 to image of Lenna 



(18.4) : Approximation 28_32 to image of Lenna 

Figure 18 
46 











(18.5) : Approximation 32_32 to image of Lenna 

Figure 18 


In the examples presented below the results obtained by applying the SWM 
are given in the frequency domain. 

Take, for instance, the approximation 8_32 to sub-image of the image of 
Lenna. The corresponding SWT can be expressed with a sequence of 64 triads. 
The third element of each triad is the coefficient considered. The first element of 
that triad is the frequency fi^x (where i = 1, 2, 3,..., 8) on the x-axis, specified 
by the first subscript of that coefficient. The second element of that triad is the 
frequency fj^y (where j = 1, 2, 3,..., 8) on the ^-axis, specified by the second 
subscript of that coefficient. The complete list of these 64 triads is as follows: 




Ci.i) 

(/2,x, 


C2,l) 


/l,2/5 

C3,l) 


/l,2/5 

C4,l) 

{f5,X') 

/by; 

C5,l) 

{f6,X'i 

/i,y; 

C'ea) 

{f7,x, 

fl,yi 

C' 7 , 1 ) 

{f8,xi 

f^,y‘'> 

^ 8 , 1 ) 


(0.5000000,0.5000000; 161.250) 
(0.5161290,0.5000000; 14.500) 
(0.5333333,0.5000000; -64.000) 
(0.5517241,0.5000000; 46.000) 
(0.5714286,0.5000000; 52.750) 
(0.5925926,0.5000000; -70.250) 
(0.6153846,0.5000000; 9.500) 
(0.6400000,0.5000000; -11.500) 
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f2,y',Ci^2) 

{f2,x, f2,y',C2,2) 
(/3,a:,/2,y;C'3,2) 
{f4,x, f2,y'i (^ 4 , 2 ) 
{f5,x, f2,y',C5^2) 
(/e,*) /2,yj ^^6,2) 
(/7,a:,/2,y;C'7,2) 
{f8,x, f2,y', (^ 8 , 2 ) 

{fl,x, f2,y\Ci^3) 

{f2,x, f3,y',C2,3) 
{f3,x, f3,y\ (^ 3 , 3 ) 
{f4,x, f3,y\ C4,3) 
(/s.a:, /s.y; (^ 5 , 3 ) 
(/c,*) f 3 ,y \ ^6,3) 
{f7,x, f3,y\C t^s) 
(/s.K, /s.y; (^ 8 , 3 ) 
(/l,*) fi,y \ ^ 1 , 4 ) 

{f2,x, fi,y\C2,4) 

{f3,x, fA,y\ (^ 3 , 4 ) 

{f4:,x, f4,y', (^ 4 , 4 ) 

{.fh^xt f^^y^ ^ 5 , 4 ) 
{f6,x, fi,y',Ce,i) 
{f7,x, fA,y\ C7,4) 

(/s.a:, fA,y\ 

{fl,x, fb,y\ Cl, 5 ) 
(/2,a:, /s.y; C2,5) 
(/s.Kj /s.y; Cs^s) 
{f^yX^ f^yy] ^4,3} 

{fbyX, fbyy\ C5,5) 

{f&yX, fbyV', Ce.s) 

{f7yX, f5yy\C7^5) 
{fSyX, fbyy\Csyb) 
{fl,x,f6yy-,Cifi) 
if2yX, /e.y; C2,6) 

{fsyx, /e.y; Cs^e) 
ifAyXt /c.yj C4^6) 

{fbyX, /e.y; C5,6) 


(0.5000000,0.5161290; -10.000) 
(0.5161290,0.5161290; -2.000) 
(0.5333333,0.5161290; 3.250) 
(0.5517241,0.5161290; -1.750) 
(0.5714286,0.5161290; -2.000) 
(0.5925926,0.5161290; 12.000) 
(0.6153846,0.5161290; -4.750) 
(0.6400000,0.5161290; -3.250) 
(0.5000000,0.5333333; -9.500) 
(0.5161290,0.5333333; 7.000) 
(0.5333333,0.5333333; -18.500) 
(0.5517241,0.5333333; 2.000) 
(0.5714286,0.5333333; 3.250) 
(0.5925926,0.5333333; -4.250) 
(0.6153846,0.5333333; 14.750) 
(0.6400000,0.5333333; 6.750) 
(0.5000000,0.5517241; -8.750) 
(0.5161290,0.5517241; -0.500) 
(0.5333333,0.5517241; 4.750) 
(0.5517241,0.5517241; -4.250) 
(0.5714286,0.5517241; -11.250) 
(0.5925926,0.5517241; 10.000) 
(0.6153846,0.5517241; -1.500) 
(0.6400000,0.5517241; 2.000) 
(0.5000000,0.5714286; -4.750) 
(0.5161290,0.5714286; -6.000) 
(0.5333333,0.5714286; -4.500) 
(0.5517241,0.5714286; -13.000) 
(0.5714286,0.5714286; -19.000) 
(0.5925926,0.5714286; 20.000) 
(0.6153846,0.5714286; -0.500) 
(0.6400000,0.5714286; -3.000) 
(0.5000000,0.5925926; -14.250) 
(0.5161290,0.5925926; 5.250) 
(0.5333333,0.5925926; 4.500) 
(0.5517241,0.5925926; -4.000) 
(0.5714286,0.5925926; 7.000) 
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{f6,X’! f6,y‘’) 

{fl.x, f6,y] 

{f8,x, fQ,y] Cs^q) 

{fl,x, fl^y] Cij) 

(/2,x, fl^y] C2j) 
{fs^x, f7,y] Csj) 
{f^.x, h.y] Ca^j) 
(/5,x, fj.y] Cs,?) 
{f6,x, f7,y] Cqj) 
{f7,x, f7,y] Cjj) 
{f8,x, f7,y] Csj) 

{fl,x, f8,y] Ci^s) 

(/2,x, f8,y] C2,8) 
{fs^x^ f8,yi C3,S) 

{fi,x, fs,y',C4^s) 
(/s.K, f8,y', C^fi) 

(/c,*) /s.y i C'e.s) 

{f7,x, fs,y', C't.s) 
(/s.Kj f8,y] Cg.s) 


(0.5925926,0.5925926; 3.500) 
(0.6153846,0.5925926; 5.250) 
(0.6400000,0.5925926; 3.750) 
(0.5000000,0.6153846; -2.250) 
(0.5161290,0.6153846; 0.000) 
(0.5333333,0.6153846; 28.750) 
(0.5517241,0.6153846; -13.750) 
(0.5714286,0.6153846; -8.750) 
(0.5925926,0.6153846; 3.500) 
(0.6153846,0.6153846; 5.000) 
(0.6400000,0.6153846; 12.500) 
(0.5000000,0.6400000; 2.750) 
(0.5161290,0.6400000; -5.250) 
(0.5333333,0.6400000; 9.750) 
(0.5517241,0.6400000; 3.500) 
(0.5714286,0.6400000; 2.250) 
(0.5925926,0.6400000; 2.250) 
(0.6153846,0.6400000; -1.750) 
(0.6400000,0.6400000; -2.250) 


The SWTs of different approximations to sub-image 19, 8 of Lenna have been 
displayed graphically in figure 19. Sub-figure 19.1 is the graphic presentation 
corresponding to the above sequence of 64 triads. (Some coefficients cannot be 
seen clearly here due to the small size of the values of their moduli.) 
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(19.1) : SWT of the approximation 8_32 to sub-image of the image of 
Lenna. The positive coefficients are shown in blue, and the negative in red. 



(19.2) : SWT of the approximation 16_32 to sub-image of the image of 
Lenna. The positive coefficients are shown in blue, and the negative in red. 

Figure 19: Graphic display of SWTs of different approximations to sub-image 
19, 8 of Lenna 
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(19.3) : SWT of the approximation 24_32 to sub-image of the image of 
Lenna. The positive coefficients are shown in blue, and the negative in red. 



(19.4) : SWT of the approximation 32_32 to sub-image of the image of 
Lenna. The positive coefficients are shown in blue, and the negative in red. 


Figure 19: Graphic display of SWTs of different approximations to sub-image 
19, 8 of Lenna 
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5 Discussion and Prospects 

Not only has it been shown how the SWM can be applied to the analysis of 1) 
functions of one variable which in the intervals where they are analyzed satisfy 
the conditions of Dirichlet; and 2) sequences of samples from different types of 
biomedical recordings (ECGs, EEGs, and EMGs); but it has also been shown 
how that method can be applied to the analysis of images. 

Eourier’s outstanding and influential contributions, along with other ap¬ 
proaches based on the development of his ideas, have played a very important 
role in a number of scientific and technological fields. In particular, mathemat¬ 
ical tools such as the Eourier series, the Eourier transform, the discrete Eourier 
transform (DET), the fast Eourier transform (EET) - an algorithm to com¬ 
pute the DET - and wavelets are used in the field of signal and image analysis. 
Perhaps in this field, the SWM can become a useful tool to complement the 
techniques provided by the Eourier approach. 

Among the advantages which the SWM offers are its simplicity and efficiency 
with which it is possible to make the necessary computations for its application, 
using software and hardware readily available; and 2) the systematic way in 
which it is applied: In all cases its use implies the solution of a system of linear 
algebraic equations which can be specified unambiguously. 

An algorithm making it possible to apply the SWM to functions of n vari¬ 
ables, where n = 1,2,3,..., will be addressed in another article. 

The authors of this paper are particularly interested in the application of 
the SWM to the analysis of signals and images essential to the field of medicine 
and expect to devote several articles to go into greater depth on the topic. 
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